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We study coherent and incoherent difFractive hadron production in high energy 
quarkonium-heavy nucleus coUisions as a probe of the gluon saturation regime of QCD. 
Taking this process as a model for pA collisions, we argue that the coherent difFractive gluon 
production, in which the target nucleus stays intact, exhibits a remarkable sensitivity to 
the energy, rapidity and atomic number dependence. The incoherent difFractive gluon pro- 
duction is less sensitive to the details of the low-x dynamics but can serve as a probe of 
fluctuations in the color glass condensate. As a quantitative measure of the nuclear effects 
on diffractive hadron production we introduce a new observable - the diffractive nuclear 
modification factor. We discuss possible signatures of gluon saturation in diffractive gluon 
production at RHIC, LHC and EIC. 

I. INTRODUCTION 

Coherent diffractive liadron production in pA collisions is a process p + A X + h + [LRG] + A, 
wliere [LRG] stands for Large Rapidity Gap. We have recently argued in (TJ [21 13] that the coherent 
diffractive hadron production can serve as a sensitive probe of the low-x dynamics of the nuclear 
matter. The coherent diffractive production exhibits a much stronger dependence on energy and 
atomic number than the corresponding inclusive process. Indeed, the diffractive amplitude is 
proportional to the square of the inelastic one. At asymptotically high energies the coherent 
diffractive events are expected to constitute a half of the total cross section, other half being 
all inelastic processes. Therefore, coherent diffraction is a powerful tool for studying the low-x 
dynamics of QCD. In particular, we advocated using the coherent diffraction as a tool for studying 
the gluon saturation [H |5] . 

The low-x region of QCD is characterized by strong gluon fields [6l|7] and can be described in the 
framework of the color glass condensate [H |9l |TOl Hi] |T2l [13]. Equations describing the color glass 
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condensate take the simplest form in the the mean-field approximation. It is valid for a color field 
of a heavy nucleus in the multicolor limit; the relevant resumation parameter is aj 

~ 1. To 

the leading order in this parameter, the nuclear color field is a coherent non-Abelian superposition 
of the color fields of single nucleons. Higher order corrections arise due to the nucleon-nucleon or 
parton-parton correlations within a nucleon. The later contribute towards the gluon saturation in 
proton. 

In all phenomenological applications of the CGC formalism, one usually relies on the mean-field 
approximations in which only the lowest order Green's functions are relevant. If rapidity interval 
Y = ln(l/x) is such that OgY < 1, then the low-x quantum evolution effects are suppressed and 
the color field can be treated quasi-classically. When OgY ~ 1, the quantum corrections become 
important and are taken into account by Balitsky-Kovchegov equation [HI [15], which is the first 
(truncated) equation in the infinite hierarchy of the coupled integro-differential equations governing 
evolution of Green's functions of all orders (JIMWLK equations P E] [TOl El HSl [l3]). 

Although corrections to the mean-field approximation, i.e. quantum fluctuations about the 
classical solution, are assumed to be small in pA collisions at RHIC their detailed phenomenological 
study is absent. In this paper we fill this gap by calculating the differential cross section for the 
incoherent diffractive gluon production which happens to be proportional to the dispersion of the 
quasi-classical scattering amplitude in the impact parameter space, as is expected in the Glauber 
theory |14] . The incoherent diffractive gluon production in pA collisions is a process 

p + A^ X + h+ [LRG] + A* , (1) 

where A* denotes the excited nucleus that subsequently decays into a system of colorless protons, 
neutrons and nuclei debris. 

Note, that the incoherent diffraction ([T]) measures fluctuations of the nuclear color field. At 
higher energies color field of proton also saturates. The mean- field approximation is not at all 
applicable in this case. A quantum fluctuations (e.g. the so-called "pomeron loops" [T6l [T71 [T8 | [T9]) 
may be the driving force of the gluon saturation in proton. This problem is of great theoretical 
and phenomenological interest. However, in spite of considerable theoretical efforts solution to this 
problem is still illusive (see |20] for a brief review of recent advances) . 

Motivated by theoretical and phenomenological interest to the incoherent diffraction as a mea- 
sure of fluctuations of the saturated nuclear color field, we derive a formula for the incoherent 
diffractive hadron production ([T| in the framework of the color glass condensate using the dipole 
model of Mueller [21] and use it for numerical study in the RHIC and LHC kinematic regions. The 
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limiting cases of large and small invariant masses of the diffractively produced system were previ- 
ously considered by many authors [22l [23l [Ml ES Ell EH Ell EH] • Inclusion of the diffractive gluon 
production was recently shown to be essential for the phenomenology of diffractive DIS off heavy 



nuclei [30]. We generalize there results for all invariant masses (satisfying (31)) and include the 
gluon evolution effects at all rapidity intervals. Unlike models based on the leading twist nuclear 
shadowing [STJ |32]) we explicitly take into account the gluon saturation effects in the nucleus. 
The paper is structured as following. In Sec. |ll] we review the Glauber approach [H] to the 



coherent and incoherent diffraction in pA collisions. In Sec. Ill we discuss generalization of this 
approached to the case of quarkonium-nucleus scattering in high energy QCD. We then turn 



in Sec. IV to the main subject of the paper ~ the diffractive hadron production - and calculate 
contributions from coherent and incoherent channels in the quasi-classical approximation neglecting 
the low-x evolution effects. The low-x evolution effects are taken into account in Sec. |Vj 



Numerical results of the diffractive cross section are presented in Sec. VI Our phenomenological 
analyses indicates that the kinematic regions of RHIC and LHC are not asymptotic as far as the 
energy dependence of the diffractive hadron production is concerned. We observed that the ratio of 
the coherent and incoherent diffractive hadron production at midrapidity in pA collisions increases 
from RHIC to LHC. This is because gluon saturation effects in proton are assumed to be small. 
Absence of this feature in experimental measurements can serve as an evidence for the onset of 
gluon saturation in proton. 

It is important to emphasize that all our results are obtained without imposing any experimen- 
tal constraints on the forward scattering angle measurements. Coherent diffraction measurements 
require much better forward scattering angle resolution than the incoherent ones. Therefore, inco- 
herent diffraction may well turn out to be the only one accessible in experiment. 

We summarize in Sec. IVIII 



II. DIFFRACTION IN THE GLAUBER MODEL 
A. Total cross section for coherent diffraction in pA 

A general approach to multiple scattering in high energy nuclear physics was suggested by 
Glauber First, consider the pp scattering and introduce the elastic pp scattering amplitude 
iTPP(s,h), where s is the center-of-mass energy squared and b is an impact parameter. At high 
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energies b is a two-dimensional vector in the transverse plane. According to the optical theorem, 

aZ = 2 J d^RerPP{s,h) . (2) 

The scattering amplitude TPP{s, b) can be written as 

rPP{s, b) = 1 - e-'^^''-''^^ , (3) 

where x^^{s, b) is a change of phase due to interaction at point b. In the following we will always 
assume the isotopic invariance of the scattering amplitudes, e.g. r^^(s,b) = rP^(s,b), where A'^ 
stands for a nucleon. 

In pA collisions momenta of nucleons can be neglected as compared to the incoming proton mo- 
mentum. Therefore, their positions b^, a = 1, . . . ,A are fixed during the interaction. If scattering 
of proton on a different nucleons is independent, then the corresponding phase shifts Xa^ {s, B — b^) 
add up: 

A 

X'^is, B, {b4) = ^ xf{s, B-ha), (4) 

a=l 

where B is proton's impact parameter. Indeed, this result holds in QCD as was demonstrated by 



A. Mueller [33], see Sec. Ill The scattering amplitude of proton p on a nucleus A reads 

A 



rP^{s,B,{ha}) = i_e-x''-^(^.B,{b4) = i_-Qg-ixr(s.B-b„) ^ i_-Q (1 _rf^(s,B -bj) . (5) 



a=l a=l 

Introduce the nucleus- averaged amplitude 



rf/(.,B) = {A,\TP^{s,B,{K})\Af), (6) 

with \Ai) being the initial and \Af) the final nucleus state. Then, by the optical theorem, the total 
pA cross section is given by 

aP^, = 2 I d'bReTPf{s,h)- (7) 

Since we neglect motion of nucleons during the interaction, the distribution of nucleons in the 
nucleus is completely specified by the thickness function r(b): 

/oo 
dzpA{h,z), (8) 
-oo 

where /9yi(b,z) is nuclear density and z is the longitudinal coordinate. The thickness function is 
often written in terms of the mean density of nucleons p = A/{^-kR^) and the nuclear profile 
function Tyi(b) as 

T{h) = pTA{h). (9) 
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It is normalized as 



j (fbT{h) = j dHpTAih) =A. 
Using these definitions we write the diagonal element of ([6| as 

rf(s,B) = J [ l[dXpTA{ha)TP^{s,B-ha) 

•' a=l 

^-jj dHarP''{s,ha)pTA{B-ha) 



(10) 



1 - e 



- / d'^ba rp^(s,b„) prA{B-b„) 



(11) 



where we assumed that 1. 

Impact parameter dependence of the pp collisions is usually parameterized as 



(12) 



where we neglected a small imaginary part of r^^(s, b) and introduced the proton profile function 
Sp{h) as 

. 9 

1 



e p . 



(13) 



Phenomenologically, Rp = \/2B ~ 1 fm, where B = 12.6 GeV ^. Since the proton radius Rp 
is much smaller than the radius of a heavy nucleus Ra, we can approximate the proton profile 



function by the delta- function in which case the overlap integral appearing in (11 1 becomes simply 
the nuclear thickness function 



/ 



d^Sp{h)pTA{B-h)^pTAiB). 



(14) 



This approximation holds when Rp <C Ra, i-e. A^^^ ^ 1. Actually, this condition follows from a 
requirement that the quasi-classical approximation holds. Indeed, in the quasi-classical approxi- 
mation as <^ 1 and a'^A^^^ ~ 1 implying that A^^^ ^ 1. 



Using (12) and (14) in (11) we derive 

rr(.,b) = i-e 



(15) 



The cross section of coherent diffraction is given by the elastic cross section 

2 f „ / 1 vN , X ^ N \ 2 



^cd(^) 



rr(^,b) 



(16) 
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B. Total cross section for incoherent diffraction in pA 



Consider incoherent diffraction of proton p on a nucleus A. In this processes the nucleus gets 
excited from the initial state \Ai) to any colorless final state \Af). It may then decay, but it is 
essential that its constituent nucleons remain intact as color objects. The corresponding cross 
section reads 



d-'BY,{Af\^'"'is,B,{ha})\A,)HAj\TP^{s,B,{K})\A, 



d''B^{Af\TP^{s,B,{K})\A,)HAf\TP^{s,B,{K})\A, 
f 

d'B\{A,\TP^{s,B,{ha})\AM^ 



{A,\ \rP^{s,B,{ha})f\A,) - \{A,\rP^{s,B,{ha})\A, 



(17) 



where in the last line we used the completeness of the set of nuclear states. We arrived at the 
well-known result that the cross section for the incoherent diffraction at a given impact parameter 
B is given by the square of the standard deviation of the scattering amplitude T'p^{s^B, {ha}) from 
its mean- field value {Ai\TP^{s,B,{ha])\A i) in a space span by impact parameters {b^}. Clearly, 
in the black-disk limit, corresponding to the asymptotically high energies s — > oo, this deviation 
vanishes because rP^(s,B,{ba}) ^ 1 for |B - b„| < Ra and is zero otherwise, see toL The 
standard deviation is a measure of quantum fluctuations near the mean-field value. 
Since, T^^ [s,{ha}) is approximately real, we have using ^ 

(rf^(5,B,{b4)' = |^i-J][i-rp^(s,B-b,)]^ 

A 

2 



A 



= 1 - 2 n [1 - r^^(5, B - ha)] + n [1 - B - ba)] 

a=l a=l 

Averaging over the nucleus and taking the large A limit we obtain 

{M\VP%,B,{ha})\^\A,) = 

1 _ 2e-^'^'^^'"^("'^-^)''^^('') + e-^'^'^[2^'"^("'^-^)-(^''^(''^"''))']''^-^(^) 



(18) 



pN f 



(19) 



where we used (12) and denoted the inelastic pN cross section as af^ [s). Substituting (19) into 



(17) and using (16) we derive 



af{s)pTA{h) 



1 - e-'^oi 



(20) 
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Elastic cross section o"^/^, which appears in (20), can be found by taking square of (12) and 
integrating over the impact parameter: 



pN 



pN\2 



87ri?2 



(21) 



where we used (13). 



It is seen from (16) and (20) that since the pN cross sections increase with energy, at asymp- 



toticahy high energies the incoherent diffraction cross section vanishes, whereas the coherent one 
reaches a half of the total cross section. This well-known conclusion is a consequence of unitarity 
of the scattering amplitude and thus is independent of interaction details. 

III. DIFFRACTION IN THE DIPOLE MODEL 



It is phenomenologically reasonable to approximate the proton light-cone wave-function (away 
from fragmentation regions) by a system of color dipoles [3l [Ml ES]. If separation of quark and 
anti-quark is small, one can apply the perturbation theory to calculate the scattering amplitude of 
quarkonium on the nucleus. It was demonstrated by Mueller in |^21j that at high energies the qqA 



forward elastic scattering amplitude takes exactly the same form as ( 15 ) with the pA cross section 



replaced by the qqA one. By virtue of the optical theorem, the total quarkonium-nucleus cross 
section reads 



arll^{s] r) = 2 / d'bNA{v, b, y) = 2 I d'b (1 - e'^''^-^ 



1 „1<lN 



{s;r)pTA(h) 



where 



iVA(r,b,y) =Rerf^(s,b;r 



(22) 



(23) 



is the imaginary part of the forward elastic qqA scattering amplitude and Y = ln(l/x) = ln(s/so) 



is rapidity with so a reference energy scale. Eq. (22 ) is called the Glauber-Mueller formula. Let us 



note that it is far from obvious that the high energy amplitude in (22) must have the same form 



as the low energy one (15). In this correspondence it is crucial that the color dipoles are identified 



as the relevant degrees if freedom at high energies. Eq. (22) holds when In(mjv-RA) <^Y <^ l/os, 
where rriAr is a nucleon mass. This condition guarantees that the coherence length Ic is much larger 



than the nuclear radius, see (30), and that the low-x gluon evolution is suppressed. In this case 



qqN cross section can be calculated in the Born approximation (two-gluon exchange) as 

^t'lf(^;r) = ^vrVxG(x,l/r2), (24) 
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with the gluon distribution function 



xG{x,l/r'): ^ 



In 



vr 



(25) 



where fj, is an infrared cutoff. We can re- write (24) in terms of the gluon saturation scale Q^q 
defined as 



m-1 



pTAih)xGix,l/r^). 



(26) 



Subscript indicates that the low-x evolution is suppressed. Eq. (25) implies that in the Born 



approximation the qqN total cross section (24) and hence the saturation scale (26) are energy 
independent. Using this definition we have 



2iv/ ^^0' 



(27) 



where Cp = {N^ - l)/(2iVc). Substituting (27) into (22) we derive another representation of the 
total quarkonium-nucleus cross section 



4 Nc ^sO 



(28) 



Integrand of (28) represents the propagator of quarkonium qq of size b through the nucleus at 



impact parameter b. Since cr^ft^ ~ and pT^(b) ~ A^^^ Eq. (28) sums up powers of a^A^^^. 
It has been shown by Kovchegov that this corresponds to the coherent scattering of qq off the 
quasi-classical field of the nucleus |7] . This is possible only if the coherence length Ic of the qq pair 
is much larger than the nuclear radius. This condition certainly holds at RHIC energies. In the 



quasi-classical approximation the coherent diffraction cross section is given by (cp. (16)) 



d^b 1 -e"3 



(29) 



Eq. (29) represents the leading contribution to the diffraction cross section in the quasi-classical 



field strength. On the contrary, the incoherent diffraction cross section, see(20), vanishes in this 



mean-field approximation since a'^^^ pT^ih) ~ a^A^^^ ~ ^ 1. Vanishing of the incoherent 



diffraction cross section can also be seen directly from (17). It is a consequence of vanishing 



relative fluctuations at large occupation numbers of classical fields. This effect has recently been 
studied within the Glauber framework (see Sec. |ll]) in [351 IMl EI] . 

Let us emphasize to avoid a possible confusion, that the semantic distinction between the 
coherent and incoherent diffraction concerns the nucleus target staying intact or breaking down 
in the collision. As far as the coherence length is concerned, at high enough energies such that 
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Ic ^ Ra the scattering is coherent for both coherent and incoherent diffraction. However, unlike 
the incoherent diffraction that can happen in the incoherent scattering at low energies (provided 
only that Ic 3> Rn), the coherent diffraction is possible only if the scattering is coherent over the 
entire nucleus. 



IV. DIFFRACTIVE GLUON PRODUCTION 
A. Coherent diffractive gluon production 

Consider coherent diffractive production of a gluon of momentum k in qqA collision. As men- 
tioned before, coherent diffraction is possible only if the coherence length Ic of the emitted gluon 
with momentum k is larger than the nucleus size Ra (in the nucleus rest frame): 

Ic='^:>Ra: (30) 

where + indicates the light-cone direction of the incoming proton. The invariant mass of the 
produced system is given by M"^ = k^/x, where x = kj^/pj^ and p is the proton momentum. 



Substituting these equations in ( 30 ) yields the following condition on the mass of the diffractive 
system: 

m2«^ = -^, (31) 
Ra RA-rrip 

where ^/s is the center-of-mass energy of the proton-nucleon collision and nip is proton mass. 

Coherent diffractive gluon production off the large nucleus was calculated in [251 EH EH] . The 
corresponding cross section reads 

X (rr^(.,b;x,y,zi)-rr(«,b;x,y) 

lf^(.,b;x,y,Z2)-rr(«,b;x,y)) , (32) 



where the qq qqG light-cone wave function ^'^'^ is given by 

^qqf \ "^'^^ /zi-X Zi-y \ /Z2-X Z2-y , , . 

$'?'?(x,y,zi,Z2) = ^ ^ - ^ • ^ - ^ . (33) 

TT^ \|zi — x|^ |zi— ypy \|Z2 — x|^ |Z2— y| 



The qq scattering amplitude Tjf^{s,h) is given by (29), while the qqG one was calculated in 



|39l SQl E] and reads 

rf^^(.,b;x,y,z.) = (^,|rf^^(.,B,{b4;x,y,z^)|^,) =l-e-5-.'r"(™.)pT'^(b) 



2^2 
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where x, y and Zo- are the transverse coordinates of quark, antiquark and gluon respectively; a = 1 
in the amphtude and o" = 2 in the complex conjugated one. For future reference, note that we can 



express the qqGN total cross section in terms of the qqN one given by (27) 

1 



qqGNf 

<ot (•5; X, y, 



2Cf 



f^tot (s; X - z^) + all^ (s; y - z^ 



c^tof (•5;x-y) 



(35) 



where we explicitly indicated the dipole size and energy dependence of the qqN total cross section. 

Depending on the relation between the gluon emission time Tq- and the interaction time in the 
amplitude and in the c.c. amplitude there are four possible products of the qq and qqG amplitudes 



appearing in the second line of ( 32 ) . Explicitly, 



{A,\T^^^{s, B, {b4; X, y) |^,)t {A,\T'i^^{s, B, {bj; x, y) |^,) 
+ (yl.in'^^^Cs, B, {b4; X, y, zi)| Ai)t B, {bj; x, y, Z2)| A 

- {A,\T'i^^{s, B, {b4; X, y) |^,)t (^.I^^G^Cs, B, {b,}; x, y , z^Mi) 

- (yl.in'^^^Cs, B, {b4; X, y , zi)| A,)t {A,\T'i^\s, B, {b,}; x, y)\A,) . 



(36) 



B. Incoherent diffractive gluon production 



To calculate the propagators in the case of incoherent diffraction we need to consider each one of 



the four cases shown in (36) and follow the by now familiar steps (17)-(20). Consider, for example. 



gluon emission before the interaction in the amplitude and in the c.c. one, i.e. < for cr = 1, 2. 
The detailed calculation is presented in Appendix. The result is (cp. (pO])) 



j;(^^|r''^"G^(.,B,{b4;x,y,zi)|A,)t(^^|rf^^(s,B,{b4;x,y,Z2)|^. 



1 qqGN , s I qqGN , ^ 1 qqGN , qqGN , s 

^2ytot («;x,y,zi)+(7tot («;x,y,Z2)- j^o-t^^t (s;x,y,zi)crt^„\ (s;x,y,Z2) 



pTA(h) 



X n - e 



(37) 



This result holds in the quasi- classical approximation 1. Therefore, contribution of elastic 



processes to (37) are of the order 



It is suppressed compared to the contribution of inelastic processes which are of order alA^/^ ~ 1. 



Note, that the contribution to the incoherent diffraction cross section given by (37) vanishes with 
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vanishing elastic qqN cross section. Thus, we expand the expression in the curly brackets of (37), 
keeping the leading elastic term. We derive 

^(^;|n5«^(s,B,{b4;x,y,zi)|A,)t(^^|n5^^(,,B,{b4;x,y,Z2)|^.) 



aZ^\s- X, y, zi) at^\s- x, y , z,) ^W^"" i^^^M P^A^) . (38) 



This approximation is valid at not too high energies satisfying 

1 / ggGAf\2 



^^rj .T,(b)«l. (39) 

Let us for simplicity consider a cylindrical nucleus for which pT^ih) = 2A/(7ri?^). Since Ra = 
A^/^Rp and taking into account that ^a^^^'^ pTA{h) = \Qlr^ for iVc > 1 (see Q and S) we 



can re- write (39) as 



Jq2^2<<2^V^ (40) 

The saturation effects in the nucleus become important when \Q'\t^ ^ 1. At rapidities Y ~ l/a^ 
the low X evolution effects give rise to the energy dependence of the saturation scale. Let Y\ be the 
rapidity at which \Q\r^ = 1 and I2 be the rapidity at which \Q\r'^ = 2A^I^ . Since Q\ tx e^'^ |l2] 
we obtain 12 — ^1 = ln(2A^/^). Using the phenomenological value A = 0.25 for the Gold nucleus 
A = 197 we find that there are Y2 — Yx ~ 6 units of rapidity between Y\ where the saturation 



effects become important and Y2 where the expansion (38) breaks down. This provides quite wide 



kinematic window in which approximation (39) holds. 

Other three cases corresponding to different relations between Tq- and can be worked out in a 
similar way. The result for the cross section of incoherent diffractive gluon production takes form 

- (^3^/.^..^..ci^...-(x,y,z.,z.)e-(-^^)pT.(b) 

^r,«ff^(s;x,y,zi)e-i'^-"(™i)''^-(^)-a,l^(s;x,y)e-i^^^^^^ 
(a,^|f^(s;x,y,Z2)e-i'^-"(™^)''^-(^)-a,^ff(s;x,y) 

(41) 



V. LOW-x EVOLUTION 



Incoherent diffraction 



The results derived in the previous section can be generalized beyond the quasi-classical level 
to include the low-x evolution. This procedure follows the general strategy developed in |43J and 
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has been used in [H |2] to derive the expression for the coherent difFractive gluon production in 
onium-nucleus cohisions. In this we section we derive an analogous expressions for the incoherent 



diffractive gluon production. For applying this general strategy it is important, that formula (38) 



can be factorized in a product of two expressions one depending on the gluon coordinate zi while 
another one on Z2. 

The low-x gluon evolution in the nucleus at large Nc is taken into account by the following 



two substitutions in (41): (i) the exponents are replaced by the forward elastic qqA scattering 



amplitude A'^^(r, b, Y) 



(42) 



Nj^{r,h,Y) evolves towards higher rapidities Y, i.e. lower x, according to the nonlinear BK evo- 



lution equation [8l |T5] from its initial condition, given by (22) at some rapidity Y = Yq; (ii) the 
factors linear in crfot^ have emerged in expansion (38 ) where the terms of higher order (i.e. multiple 



scattering) in the elastic amplitude were neglected. Therefore, they are replaced by the forward 
elastic qqN scattering amplitude Np{r, b, Y) that evolves according to the linear BFKL equation 



It is defined similarly to (23) (see (112 



iVp(r,b,y) =Rer'?*(s,b;r) = -^^^(s; r) 5p(b) 



(43) 



Recall, that in the heavy nucleus environment the impact parameter dependence of qqp cross section 



can be neglected, see (14). Thus, using (13) we have 



a?|f(s;r)^27ri?^iVp(r,0,y) 



(44) 



in agreement with the optical theorem. The factors linear in the qqGN cross section are replaced 



similarly using (35). Substituting (42) and (44) into (41) using (35) (in the multicolor Nc ^ 1 
limit) yields 



j d%(fzi <fz2 y, Zi, Z2) pTA{h) e-^k.(zx-Z2) 



dcrm{k,y) _ 1 TrRp 
d?kdy ~ (2^)2 2 

X I [1 - iV^(zi - X, b, y)] [1 - Na{2i - y , b, y)] [Afp(zi - x, 0, y) + iVp(zi - y , 0, y)] 



[l_Ar^(x-y,b,y)]iVp(x-y,0,y) 



X I [1 - Na{22 - X, b, y)] [1 - Na{22 - y , b, y)] [A^p(z2 - x, 0, y) + iVp(z2 - y , 0, y)] 



[l_Ar^(x-y,b,y)]iVp(x-y,0,y) 



(45) 
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We can write (45) in a more compact form introducing the following two-dimensional vector 
Iid(x - y, b, y; k) = / <fz f - ^^"i e-**^'^ 

X I [1 - iV^(z - X, b, y)] [1 - iV^(z - y , b, y)] [A/p(z - x, 0, y) + iVp(z - y, 0, y)] 



[1 - 7Va(x - y , b, y)] iVp(x - y , 0, y) | . 



(46) 



Then, 



y d26prA(b)|IiD(x-y,b,y;k)|2. (47) 



d'^kdy 7r2 (27r)2 2 

So far we have discussed the gluon evolution in the rapidity interval y between the emitted gluon 
and the target nucleus. The low-x evolution also occurs in the rapidity interval Y — y between 
the incident quarkonium qq and the emitted gluon (we now denote by Y the rapidity between 
the quarkonium and the nucleus and by y the inclusive gluon rapidity). The low-x evolution 



in the incident quarkonium is taken into account by convoluting (47) with the dipole density 
ni(r, r', B — h,Y — y), where r = x — y: 



d'^kdy 7r2 2(27 



^ Jd%J d^B j (iVni(r,r',B-b,y-y)prA(b)|IiD(r,b,y;k)|2. 

(48) 

The dipole density ni(r,r,B — b,y — y) has the meaning of the number of dipoles of size r' at 
rapidity Y — y and impact parameter b generated by evolution from the original dipole r having 
rapidity Y and impact parameter B [211 . It obeys the BFKL equation |44| 35]. As the interval 
Y — y increases the quarkonium wave function involves increasing number of dipoles. Although the 
dipole size distribution diffuses according to the BFKL equation, the typical dipole size is much 
smaller than the nuclear radius Ra, i.e. |B — b| ^ h. Integrating over B — b in this approximation 



we derive 

d(JiT){k,y) asCp T^^lp 



d^kdy 7r2 2(27r)2 

where we defined [U [2] 



j dH I d\'np{r,r',Y-y)pTA{h)\hG{r,h,y;k)\\ (49) 



np(r, r', y) = J ni{r, r', b', Y) dH' . (50) 



Eqs. (49) and (46) constitute the main result of this paper. 
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B. Coherent diffraction 

The cross section for the coherent difFractive gluon production including the low-j; evolution 
was derived in [H |2] . It can be written similarly to ( |49| ) [U [2] : 

dacY){k,y) _ UsCf 1 [ f , 
d?k dy 

where 



^(2^ / '^^^ I d'r'np{r,r',Y -y)\IcD{r,h,y;k)\\ (51) 



Icd(x - y , b, y;k) = I d'z ( - ) e~^^-^ 



z — |z — y I 

X <{ - Na{z - X, b, y) - Na{z - y , b, y) + Na{x - y , b, y) 



NA{z-^,h,y)NA{z-y,h,y)}. (52) 



In [H |2] we presented a detailed analytical and numerical analysis of the the coherent difFractive 
gluon production. 

C. Logarithmic approximation 

To obtain a working model for numerical calculations, it is useful to estimate IiD(r,b,y;k) in 
the logarithmic approximation. First, let us change the integration variable w = z — y and define 
an auxiliary function [H |2] 

QiD(r', y; k) = e^^^^ | [l - TV^Cr', b, y)] Np{r' , 0, y) 

- [l-iVA(w-r',b,y)] [l-iVA(w,b,y)] [iVp(w - r', 0, y) + 7Vp(w, 0, y)] | (53) 



With this definition (46) becomes 



IiD(r, b, y; k) = -e-''^ '' iVkQm{r', y; k) + e^'^''' iVkQMr', y, k) . (54) 
In the logarithmic approximation Qij){r' ,y;k.) is a real function^j ,2^. Therefore, we can write 

|IlD(r,b,y;k)|2 = 4sin2 (^k • r') (VkQiD(r', y; k))^ . (55) 



Notice that the integral over z in (46) is dominated by the region w < 1/k since otherwise the 
integrand is a rapidly oscillation function. Restricting integration to this region yields results valid 
as long as the scales 1/r', Qg and k are strongly ordered [B 12]. Operating with Vk on QiD(r', y; k) 
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yields 

VkQiD(r',y;k) = -i^^QiD(r',y;k) = -2vrp| [l - NA{r',h,y)] Np{r',0,y) 
- [1 - NA{k k-^ - r', b, y)] [1 - NA{k k-\ b, y)] [iVp(k - r', 0, y) + Nj,{\<i k-\ 0, y)] j^SG) 
Analogously, in the coherent diffraction case we define 

QcD(r',y;k) = j ^e*k.w|^^(^_^/^b,y)+iV^(w,b,2/)-7V^(r',b,y) 

-iVA(w-r',b,y)ArA(w,b,y)|. (57) 
Then, in the logarithmic approximation, 

VkQcD(r',2/;k) = -2^^| -Af^(r',b,y) + iVA(kA:-2-r',b,y)+iV^(kfc-2,b,y) 



A:2 ^ 

iVA(k A:-2 - r', b, y) Na{\^ k~\ h,y)\ . (58) 



Similarly (55) becomes 



|IcD(r, b, y- k)|2 = 4sin2 {\ k • r') (VkQcD(r', y; k))^ . (59) 



VI. NUMERICAL CALCULATIONS 

In our previous paper ^ we justified using the quarkonium-nucleus qqA scattering as a model of 
pA collisions. It also serves as the QCD ingredient for the deep inelastic eA scattering. Therefore, 
we extend the use of this model in this section for the incoherent diffraction. 

We fix the quarkonium size at r = 0.2 fm. Dependence of all cross sections on r has been 
discussed at length in [H |2l [3] where the interested reader is readily referred. Since our main 
goal in this section is to illustrate the key features of the various diffraction channels, rather than 
giving a detailed quantitative analyses, we approximate the nuclear profile by the step function 
(cylindrical nucleus). 

It is clear from the discussion in the previous Section that in order to calculate the diffractive 
cross sections we need only to specify the dipole-nucleus scattering amplitude A^yi(r, b, y) and the 
dipole density np{r,r',y). The dipole-proton scattering amplitude Np{r,h,y) is obtained from 
Na^ti^, b, y) in the limit A — > 1. 

For the forward elastic dipole-nucleus scattering amplitude A'^^(r, b,y) we use the KKT model 
|46| . which is based on analytical analysis in [47J. It successfully describes the inclusive hadron 
production in the RHIC kinematic region (at the central and forward rapidities). 
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For the dipole density we use the leading BFKL solution in the diffusion approximation: 



np{r, r',Y - y) 



1 1 



vr 



,(OLp-l){Y-y) ^ lAC(3)asd(Y-y) _ 



(60) 



27r2 rr'Y 14C(3)qs d - y) 
Parameter d is equal to unity in the LO BFKL. To obtain the hadron diffractive cross section we 



convoluted the gluon cross sections (49) and (51) with the LO pion fragmentation function given 
in HHl. 



A convenient way to study the nuclear dependence of particle production is to consider the 
nuclear modification factor defined as follows 



d^Ttrdy 



(61) 



We calculate the diffractive gluon production in pp collisions, which is required as a baseline 



for the calculation of the nuclear modification factor (61), by setting A = 1 m the formula for 



the corresponding cross section in pA collisions. The nuclear modification factor is defined in 
such a way that a completely incoherent scattering would yield i?j^(A;2^,y) = 1. The results of 
the calculations are exhibited in Fig. [T] We see that the nuclear modification factor is rapidly 
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FIG. 1: (Color online). Nuclear modification factor for (a) coherent diffraction, (b) incoherent diffraction 
for RHIC at different rapidities. 



decreasing with rapidity for the coherent diffraction case and is almost independent of rapidity for 
the incoherent case. This trend continues even at the forward LHC energies In contrast to 
inclusive hadron production that saturates already in the forward rapidities at RHIC |46| H9] , we 
expect that the coherent diffraction cross section saturates at much higher energies and rapidities, 
perhaps at the forward rapidities at LHC (see [3i for an extensive discussion). As Fig. [T] implies, 
the incoherent diffractive gluon production is saturated already in the central rapidity region at 
RHIC. 
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We observe in Fig. [l] a parametric enhancement by a factor ~ A^^^ of the coherent diffractive 
cross section with respect to the incoherent one. This is a benchmark of the classical gluon field 
of the nucleus in which quantum fluctuations are suppressed by A^^^. This feature is also seen in 
Fig. |2] where we show the ratio of the cross sections for coherent and incoherent diffractive gluon 
production for pA collisions at different energies. To illustrate the atomic number dependence of 
this ratio we show three cases: A = 200, A = 100 and A = 20. 
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FIG. 2: (Color online). Ratio of the cross sections for coherent and incoherent diffractive gluon production 
in pA collisions at midrapidity: (a) RHIC and (b) LHC. Solid blue line: A = 200, dashed red line: A = 100 
and dotted green line: A = 20. 



In Fig. |2] we can see that the fraction of the incoherent diffractive events increases from RHIC 
to LHC. This happens because the gluon saturation effects in proton are still not strong enough 
to unitarize the cross section. Remember, this allowed us to expand the scattering amplitude as 



shown in (38). Our estimate (40) shows that this approximation is valid for about six units of 



rapidity starting from the rapidity yi ~ 1 at RHIC. Therefore, our conclusion that the fraction 
of incoherent diffractive events increases holds up to the rapidity 2/2 — 3 at the LHC. As soon 
as the gluon distribution of the incident proton saturates, the incoherent diffraction cross section 



vanishes, as can be seen in (20). This effect however is not included in our present model. 



Finally, it may be challenging to experimentally distinguish the coherent and incoherent diffrac- 
tive cross sections. In part it is related to the difficulty of performing measurements at very forward 
angles. One usually measures an event with large rapidity gap between the produced hadronic sys- 
tem and the target remnants in the forward direction (that may or may not be an intact nucleus). 
In such a likely case we consider the nuclear modification factor for a sum of coherent and in- 
coherent diffractive channels. The result is displayed in Fig. |3j We can see that the predictions 
for the total diffractive cross section are similar to those of the inclusive hadron production 
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Therefore, it is very important to experimentally separate contributions of coherent and incoher- 
ent diffraction. The former is especially useful for studying the low x QCD due to its strong and 
non-trivial energy, rapidity and atomic number dependence. 
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FIG. 3: (Color online). Nuclear modification factor for the sum of coherent and incoherent diffractive cross 
sections for (a) RHIC, (b) LHC at different rapidities. 



Let us also mention that electromagnetic interactions also contribute to diffraction processes in 
pA collisions. If they become large, as argued in [50J, it may present an additional experimental 
challenge. We intend to analyze contribution of the e.m. interactions in a separate publication. 

To conclude this section, we would like to emphasize that Figs. 1|2|3 are an illustration of 
general features expected in various diffractive channels. As we have mentioned more than once 
in this paper, application of realistic experimental cuts, which are unique for every experiment, 
will significantly change the absolute values of the cross sections and relative importance of the 
coherent and the incoherent channels. A dedicated study is required in each case, although the 
key features displayed in Figs. 1|2|[3 will perhaps remain unchanged. 



VII. SUMMARY 



To summarize, we calculated the cross section for the incoherent diffractive gluon production 
in qqA collisions, which is a prototype of pA collisions at RHIC and LHC and "f* A ones at EIC. 
We took into account the gluon saturation, i.e. color glass condensate, effects. We then compared 
prediction of diffractive hadron production in coherent (calculated previously in [H 121 [3]) and 
incoherent diffraction. 

Coherent diffractive gluon production has a very characteristic energy, rapidity and atomic 
number dependence which makes it a powerful tool for studying the leading mean-field contribution 
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to the color glass condensate. Incoherent diffraction arises from contributions to the gluon field 
correlations beyond the mean-field approximation. Its experimental study can reveal the dynamics 
of quantum fluctuations in the CGC. Unlike the nuclear modification factor for coherent diffractive 
gluon production the nuclear modification factor for incoherent diffraction is not expected to exhibit 
a significant rapidity and energy dependence. Ratio of the coherent and incoherent inclusive 
diffractive cross sections is predicted to increase from RHIC to LHC if the gluon saturation effects 
in proton are small. Otherwise, the ratio will decrease approaching the unitarity limit (i.e. zero). 

Finally, one may consider measuring a diffractive event with large rapidity gap but without 
distinguishing contributions of coherent and incoherent components. Although in this case many 
interesting features of coherent and incoherent channels get averaged out, the corresponding cross 
section as well as the nuclear modification factor are expected to display a non-trivial behavior 
as a function of energy, atomic number and transverse momentum. This behavior is sensitive to 
the underlying parton dynamics and thus can serve as a discriminator of different models. We 
are certain that studying diffraction at RHIC, LHC and EIC can become an important tool in 
accessing the detailed structure of QCD at low-x. 
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APPENDIX A 



In this Appendix we derive Eq. (37). For notation simplicity we replace the coordinates x, y, 



in the scattering amplitudes by subscript a, with a = 1,2. We have 

j;(A;|rf ^^(., B, {b4)| A,)t(Ay |rf ^^(., B, {b4)| A,) = 
(yi,|rf ^^(s, B, {b4) rf ^^(s, B, {b4)| A,) 

-(^,|rf^^(s,B,{b4)|A,)(^,|rf^^(s,B,{b4)|^,), (ai) 
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Similarly to (18) we write 



rf^^(.,B,{bj)rf^^(.,B,{b4) = i-fj [i-rf^^(s,B-b,) 



a=l 



n [i - rf «^(., B - b,)] + n [i - rf'^'^is, B - b„)] [i - rf b - k)] (A2) 



a=l 



a=l 



Averaging over the nucleus 



(A,|rf ^^(5, B, {b4) rf ^^(s, B, {K})\A,) 

= 1 _ g-/rf2brf«^(s,B-b)pT^(fe) _ g-/d2brf«^(s,B-b)pT^(fe) 

/rf2b[rf^^{s3-b)+rf'^'^(s,B-b)-rf'^^(s,B-b)rf'^'^(s,B--b)]pTA(b) 



(A3) 



and subtracting the coherent part 



1 - rf ^^(., B, {b4) ) ( 1 - Tf^^is, B, {b4) 



■^qqGA, 



(A4) 



we obtain 



j;(^^|rf«^(.,B,{b4)|^,)t(^^|rf^^(.,B,{b4)|^,) = 

g- / d^fc [rf °^(s,B-b)+rf'^'^{s,B-b)-r?'^^(s,B-b)rf '''^{s,B-b)]prA(6) 

/ '^'^ [r?'?'^'^(s,B-b)+rf °^(s,B-b)]pT^(6) 

^ J d2ft [rf °^(s,B-b)+r«'?''^{s,B-b)-r«''^'^(s,B-b)r«'?'3^(s,B-b)]pTA(b) 



X n - e- 



-/d2brf''^{s,B-b)rf°^(s,B-b)pTA(b)\ 



4^ t°t 



'(.,l)a«„f^(s,2) 



X n -e 



(A5) 
(A6) 

(A7) 



as advertised. In the last line we employed the approximation Rp <^ Ra which allows to write 



using an analogue of (12| 



d^rf^^is, B - b)rf {s, B - b) pTA{b) 



qqGN I 



,(b-B)2 



qqGA, 



pTa{B) I (fh'e 



\'^!!^%,l)<rl!^%,2)^^pTA(B) 



(A8) 
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where b' = b - B and b'^ < b^, B^. 
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